Genetic lineage tracing identifies adaptive mechanisms of pancreatic islet ß cells in various mouse models of diabetes with distinct age of initiation.

During the pathogenesis of type 1 diabetes (T1D) and type 2 diabetes (T2D), pancreatic islets, especially the ß cells, face significant challenges. These insulin-producing cells adopt a regeneration strategy to compensate for the shortage of insulin, but the exact mechanism needs to be defined. High-fat diet (HFD) and streptozotocin (STZ) treatment are well-established models to study islet damage in T2D and T1D respectively. Therefore, we applied these two diabetic mouse models, triggered at different ages, to pursue the cell fate transition of islet ß cells. Cre-LoxP systems were used to generate islet cell type-specific (α, ß, or δ) green fluorescent protein (GFP)-labeled mice for genetic lineage tracing, thereinto ß-cell GFP-labeled mice were tamoxifen induced. Single-cell RNA sequencing (scRNA-seq) was used to investigate the evolutionary trajectories and molecular mechanisms of the GFP-labeled ß cells in STZ-treated mice. STZ-induced diabetes caused extensive dedifferentiation of ß cells and some of which transdifferentiated into a or δ cells in both youth- and adulthood-initiated mice while this phenomenon was barely observed in HFD models. ß cells in HFD mice were expanded via self-replication rather than via transdifferentiation from α or δ cells, in contrast, α or δ cells were induced to transdifferentiate into ß cells in STZ-treated mice (both youth- and adulthood-initiated). In addition to the re-dedifferentiation of ß cells, it is also highly likely that these "α or δ" cells transdifferentiated from pre-existing ß cells could also re-trans-differentiate into insulin-producing ß cells and be beneficial to islet recovery. The analysis of ScRNA-seq revealed that several pathways including mitochondrial function, chromatin modification, and remodeling are crucial in the dynamic transition of ß cells. Our findings shed light on how islet ß cells overcome the deficit of insulin and the molecular mechanism of islet recovery in T1D and T2D pathogenesis.

Interaction between the GCKR rs1260326 variant and serum HDL cholesterol contributes to HOMA-ß and ISIMatusda in the middle-aged T2D individuals.

This study aims to investigate the correlations between islet function/ insulin resistance and serum lipid levels, as well as to assess whether the strength of such correlations is affected by the GCKR rs1260326 variant in healthy and T2D individuals. We performed an oral glucose tolerance test (OGTT) on 4889 middle-aged adults, including 3135 healthy and 1754 T2D individuals from the REACTION population study in the Nanjing region. We also measured their serum lipid levels and genotyped for rs1260326. We found that serum high-density lipoprotein (HDL) cholesterol and triglyceride (TG) levels were independently correlated with indexes of islet function (HOMA-ß and IGI [insulinogenic index]) and insulin resistance (HOMO-IR and ISIMatsuda) in both healthy and T2D individuals. The correlations were significantly decreased in T2D individuals, with significant heterogeneities compared to healthy controls (I2 > 75%, Phet < 0.05). Although no correlation was observed between serum total cholesterol (TC) level and islet function/ insulin resistance in healthy controls, significant correlations were found in T2D individuals, with significant heterogeneity to healthy controls in the correlation with ISIMatsuda(I2 = 85.3%, Phet = 0.009). Furthermore, we found significant interactions of the GCKR rs1260326 variant for the correlations between serum HDL cholesterol and HOMA-ß/ISIMatsuda in T2D subjects (P = 0.015 and 0.038, respectively). These findings illustrate that distinct correlations between serum lipid levels and islet function/ insulin resistance occurred in T2D subjects compared to healthy individuals. Common gene variants, such as rs1260326, might interact substantially when studied in specific populations, especially T2D disease status.

Single-cell RNA sequencing combined with single-cell proteomics identifies the metabolic adaptation of islet cell subpopulations to high-fat diet in mice.

AIMS/HYPOTHESIS: Islets have complex heterogeneity and subpopulations. Cell surface markers representing alpha, beta and delta cell subpopulations are urgently needed for investigations to explore the compositional changes of each subpopulation in obesity progress and diabetes onset, and the adaptation mechanism of islet metabolism induced by a high-fat diet (HFD). METHODS: Single-cell RNA sequencing (scRNA-seq) was applied to identify alpha, beta and delta cell subpopulation markers in an HFD-induced mouse model of glucose intolerance. Flow cytometry and immunostaining were used to sort and assess the proportion of each subpopulation. Single-cell proteomics was performed on sorted cells, and the functional status of each alpha, beta and delta cell subpopulation in glucose intolerance was deeply elucidated based on protein expression. RESULTS: A total of 33,999 cells were analysed by scRNA-seq and clustered into eight populations, including alpha, beta and delta cells. For alpha cells, scRNA-seq revealed that the Ace2low subpopulation had downregulated expression of genes related to alpha cell function and upregulated expression of genes associated with beta cell characteristics in comparison with the Ace2high subpopulation. The impaired function and increased fragility of ACE2low alpha cells exposure to HFD was further suggested by single-cell proteomics. As for beta cells, the CD81high subpopulation may indicate an immature signature of beta cells compared with the CD81low subpopulation, which had robust function. We also found differential expression of Slc2a2 in delta cells and a potentially stronger cellular function and metabolism in GLUT2low delta cells than GLUT2high delta cells. Moreover, an increased proportion of ACE2low alpha cells and CD81low beta cells, with a constant proportion of GLUT2low delta cells, were observed in HFD-induced glucose intolerance. CONCLUSIONS/INTERPRETATION: We identified ACE2, CD81 and GLUT2 as surface markers to distinguish, respectively, alpha, beta and delta cell subpopulations with heterogeneous maturation and function. The changes in the proportion and functional status of islet endocrine subpopulations reflect the metabolic adaptation of islets to high-fat stress, which weakened the function of alpha cells and enhanced the function of beta and delta cells to bring about glycaemic homeostasis. Our findings provide a fundamental resource for exploring the mechanisms maintaining each islet endocrine subpopulation's fate and function in health and disease. DATA AVAILABILITY: The scRNA-seq analysis datasets from the current study are available in the Gene Expression Omnibus (GEO) repository under the accession number GSE203376.

Rs864745 in JAZF1, an Islet Function Associated Variant, Correlates With Plasma Lipid Levels in Both Type 1 and Type 2 Diabetes Status, but Not Healthy Subjects.

Objective: This study aims to reveal the association between JAZF1 rs864745 A>G variant and type 2 diabetes (T2D), type 1 diabetes (T1D) risk, and their correlation with clinical features, including islet function, islet autoimmunity, and plasma lipid levels. Methods: We included 2505 healthy controls based on oral glucose tolerance test (OGTT), 1736 unrelated T2D, and 1003 unrelated autoantibody-positive T1D individuals. Binary logistic regression was performed to evaluate the relationships between rs864745 in JAZF1 and T2D, T1D, and islet-specific autoantibody status under the additive model, while multiple linear regression was used to assess its effect on glycemic-related quantitative traits and plasma lipid levels. Results: We did not find any association between rs864745 in JAZF1 and T2D, T1D, or their subgroups (All P > 0.05). For glycemic traits, we found that the G allele of this variant was significantly associated with higher 120 min insulin level, insulinogenic index (IGI), corrected insulin response (CIR), and acute insulin response (BIGTT-AIR) (P = 0.033, 0.006, 0.009, and 0.016, respectively) in healthy individuals. Similar associations were observed in newly diagnosed T2D but not T1D individuals. Although this variant had no impact on islet autoimmunity (All P > 0.05), significant associations with plasma total cholesterol (TC) and low-density lipoprotein (LDL) level stratified by JAZF1 rs864745 variant were observed in the disease status of T2D (P = 0.002 and 0.003) and T1D (P = 0.024 and 0.009), with significant heterogeneity to healthy individuals. Conclusions: The common JAZF1 rs864745 variant contributes to islet function and lipid metabolism, which might be put into genetic risk scores to assess the risk of related clinical features.

Genome-Wide Identification of N6-Methyladenosine Associated SNPs as Potential Functional Variants for Type 1 Diabetes.

Objectives: N6-methyladenosine (m6A) is essential in the regulation of the immune system, but the role that its single nucleotide polymorphisms (SNPs) play in the pathogenesis of type 1 diabetes (T1D) remains unknown. This study demonstrated the association between genetic variants in m6A regulators and T1D risk based on a case-control study in a Chinese population. Methods: The tagging SNPs in m6A regulators were genotyped in 1005 autoantibody-positive patients with T1D and 1257 controls using the Illumina Human OmniZhongHua-8 platform. Islet-specific autoantibodies were examined by radioimmunoprecipitation in all the patients. The mixed-meal glucose tolerance test was performed on 355 newly diagnosed patients to evaluate their residual islet function. The functional annotations for the identified SNPs were performed in silico. Using 102 samples from a whole-genome expression microarray, key signaling pathways associated with m6A regulators in T1D were comprehendingly evaluated. Results: Under the additive model, we observed three tag SNPs in the noncoding region of the PRRC2A (rs2260051, rs3130623) and YTHDC2 (rs1862315) gene are associated with T1D risk. Although no association was found between these SNPs and islet function, patients carrying risk variants had a higher positive rate for ZnT8A, GADA, and IA-2A. Further analyses showed that rs2260051[T] was associated with increased expression of PRRC2A mRNA (P = 7.0E-13), and PRRC2A mRNA was significantly higher in peripheral blood mononuclear cell samples from patients with T1D compared to normal samples (P = 0.022). Enrichment analyses indicated that increased PRRC2A expression engages in the most significant hallmarks of cytokine-cytokine receptor interaction, cell adhesion and chemotaxis, and neurotransmitter regulation pathways. The potential role of increased PRRC2A in disrupting immune homeostasis is through the PI3K/AKT pathway and neuro-immune interactions. Conclusion: This study found intronic variants in PRRC2A and YTHDC2 associated with T1D risk in a Chinese Han population. PRRC2A rs2260051[T] may be implicated in unbalanced immune homeostasis by affecting the expression of PRRC2A mRNA. These findings enriched our understanding of m6A regulators and their intronic SNPs that underlie the pathogenesis of T1D.

Phosphoproteome reveals molecular mechanisms of aberrant rhythm in neurotransmitter-mediated islet hormone secretion in diabetic mice.

BACKGROUND: Acetylcholine (ACh) and norepinephrine (NE) are representative neurotransmitters of parasympathetic and sympathetic nerves, respectively, that antagonize each other to coregulate internal body functions. This also includes the control of different kinds of hormone secretion from pancreatic islets. However, the molecular mechanisms have not been fully elucidated, and whether innervation in islets is abnormal in diabetes mellitus also remains unclear. METHODS AND RESULTS: Immunofluorescence colocalization and islet perfusion were performed and the results demonstrated that ACh/NE and their receptors were highly expressed in islet and rapidly regulated different hormones secretion. Phosphorylation is considered an important posttranslational modification in islet innervation and it was identified by quantitative proteomic and phosphoproteomic analyses in this study. The phosphorylated islet proteins were found involved in many biological and pathological processes, such as synaptic signalling transduction, calcium channel opening and insulin signalling pathway. Then, the kinases were predicted by motif analysis and further screened and verified by kinase-specific siRNAs in different islet cell lines (αTC1-6, Min6 and TGP52). After functional verification, Ksr2 and Pkacb were considered the key kinases of ACh and NE in insulin secretion, and Cadps, Mlxipl and Pdcd4 were the substrates of these kinases measured by immunofluorescence co-staining. Then, the decreased expression of receptors, kinases and substrates of ACh and NE were found in diabetic mice and the aberrant rhythm in insulin secretion could be improved by combined interventions on key receptors (M3 (pilocarpine) or α2a (guanfacine)) and kinases (Ksr2 or Pkacb). CONCLUSIONS: Abnormal innervation was closely associated with the degree of islet dysfunction in diabetic mice and the aberrant rhythm in insulin secretion could be ameliorated significantly after intervention with key receptors and kinases in the early stage of diabetes mellitus, which may provide a promising therapeutic strategy for diabetes mellitus in the future.

Interferon-α promotes MHC I antigen presentation of islet ß cells through STAT1-IRF7 pathway in type 1 diabetes.

Type 1 diabetes (T1D) is a T cell-mediated autoimmune disease. Increased incidence of T1D was reported in patients receiving IFN-α treatment. However, the exact mechanisms of IFN-α that facilitate the pathogenesis of T1D are not fully understood. To explore the mechanism of IFN-α on the immune system and islets, non-obese diabetic (NOD) mice were injected with IFN-α and the progression of autoimmune insulitis was assessed by haematoxylin and eosin (HE) staining, immunohistochemical and flow cytometry analysis. Transcriptional profiling of islets treated with IFN-α was explored by RNA-seq. IFN-α induced antigen presentation was evaluated by qRT-PCR, western blot and immunofluorescence, and key transcription factors were inhibited by small interfering RNAs (siRNAs). Our data show that IFN-α contributed to the progression of autoimmune insulitis in NOD mice by promoting the proliferation of CD8+ T cells. IFN-α upregulated antigen presentation related genes MHC I, TAP1, B2M, PSMB8, NLRC5 and transcriptional regulator STAT1, STAT2, IRF7 at a time and dose-dependent manner. The silence of STAT1 or STAT2 both weakened IFN-α-induced increase of antigen presenting related molecules. IRF7 was also merely influenced by STAT1 silence. The knockdown of IRF7 decreased the IFN-α induced expressions of TAP1, PSMB8 and MHC I and prevented the expression of STAT2 but not STAT1. Our study demonstrated that STAT1-IRF7-MHC I complex axis were crucial for IFN-α signalling in islets, and created positive feedback through IRF7-STAT2 cascade amplifying signals which accelerated the process of T1D.

Dynamic Number and Function of IL-10-Producing Regulatory B Cells in the Immune Microenvironment at Distinct Stages of Type 1 Diabetes.

The critical role of IL-10-producing B cells (B10 cells) with a unique CD1dhiCD5+ phenotype in suppressing autoimmune responses and relieving inflammation has been demonstrated in several models of autoimmune diseases. However, the regulatory role of B10 cells in T cell-mediated autoimmune responses during the natural history of type 1 diabetes is unclear. In this study, we used the NOD mouse model of autoimmune diabetes to clarify the changes and potential mechanisms of B10 cells for disease. Compared with B10 cells present in the 4-wk-old normoglycemic NOD mice, the frequency of B10 cells was increased in the insulitis and diabetic NOD mice, with the highest proportion in the insulitis NOD mice. The changes in the relative number of B10 cells were most pronounced in the pancreas-draining lymph nodes. The pathogenic T cells, including Th1 and Th17 cells, remarkably increased. The assays in vitro showed that B10 cells in the NOD mice did not inhibit the proliferation of CD4+CD25- T cells. They also had no regulatory effect on IFN-γ and IL-4 secretion or on Foxp3 expression of T cells. B10 cells suppressed T cell-mediated autoimmune responses via an IL-10-dependent pathway. In contrast, B10 cells in the NOD mice exhibited a significant reduction in IL-10 production. In summary, a defect in the number and function of B10 cells may participate in the development and progression of type 1 diabetes.

Constructing a metabolic integral score model for the quantification of metabolic dysfunction and tendency.

BACKGROUND AND AIMS: The binary nature of metabolic syndrome (MetS) cannot quantitatively describe the severity of metabolic abnormalities. We aim to establish a metabolic integral score (MIS) model to quantify the severity and polarity of metabolic disorders and their relationship with insulin sensitivity and secretion. METHODS AND RESULTS: We performed factor analysis on 9950 participants from a cross-sectional study conducted in China. The MIS model was established using 10 variables including body mass index (BMI), waist circumference, hip circumference, glycosylated hemoglobin (HbA1c), fasting and 2-h plasma glucose (FPG, 2h-PG), systolic blood pressure (SBP), diastolic blood pressure (DBP), high-density lipoprotein (HDL) and triglyceride (TG) levels. Four common factors were identified as "glucose factor," "obesity factor," "blood pressure factor," and "lipid factor," respectively, in MIS model (KMO = 0.755, P < 0.001). MIS = 0.433 × Factor 1 + 0.267 × Factor 2 + 0.172 × Factor 3 + 0.128 × Factor 4. Insulin sensitivity and ß-cell function decreased with the increase of MIS (P < 0.001). We classified four metabolic tendencies according to factor quartiles. Individuals in Tendency 1 (severe hyperglycemia) had the worst ß-cell function. Tendency 3 (severe hypertension) had the best insulin sensitivity. Tendency 4 (severe dyslipidemia) had preferable ß-cell function (P < 0.05). CONCLUSIONS: Our MIS model provides a quantitative scoring system to assess various patterns of metabolic abnormality that indicate different underlying pathophysiology.

Differences of Circulating CD25hi Bregs and Their Correlations with CD4 Effector and Regulatory T Cells in Autoantibody-Positive T1D Compared with Age-Matched Healthy Individuals.

Circulating CD25hi B cells, a subset of regulatory B cells in humans, are closely related to inflammation and autoimmune diseases. This study is aimed at investigating the alternation of CD25hi Bregs and their correlation with CD4 effector and regulatory T cells in T1D individuals. We included 68 autoantibody-positive T1D and 68 age-matched healthy individuals with peripheral blood mononuclear cells (PBMCs) and assessed them with CD25hi Bregs and CD4 effector or regulatory T cells by flow cytometry. Here, we demonstrate that the frequency of CD25hi Bregs was significantly decreased in T1D subjects (P = 0.0016), but they were not affected by disease status (age at T1D diagnosis or duration) or T1D risk loci (rs2104286 or rs12251307) in IL2RA (all P > 0.05). Moreover, higher IgD (P = 0.043) and lower CD27 (P = 0.0003) expression was observed in CD25hi Bregs of T1D individuals, but not the expression of IgM, CD24, or CD38 (all P > 0.05). Although there was no correlation between CD25hi Bregs and CD4 effector T cell subsets in either T1D or healthy individuals (all P > 0.05), we found a positive correlation between CD25hi Bregs and CD4 Tregs in healthy controls (Sp. r = 0.3544, P = 0.0249), which disappeared in T1D subjects (Sp. r = 0.137, P = 0.401). In conclusion, our results suggest that decreased CD25hi Bregs and alternation of their phenotypes are features of T1D regardless of disease duration and T1D genetic risk loci, and an impaired balance between CD25hi Bregs and CD4 Tregs might contribute to the pathogenesis of T1D.

The common rs13266634 C > T variant in SLC30A8 contributes to the heterogeneity of phenotype and clinical features of both type 1 and type 2 diabetic subtypes.

AIMS: T2D and T1D are phenotypically heterogeneous. This study aims to reveal the relationship between the common SLC30A8 rs13266634 variant and subgroups of T2D and T1D and their clinical characteristics. METHODS: We included 3158 OGTT-based healthy controls, unrelated 1754 T2D, and 1675 autoantibody-positive T1D individuals. The associations between rs13266634 and subtypes of T2D, T1D, autoantibody status and glycemic-related quantitative traits were performed by binary logistic regression analysis under the additive model and multiple linear regression with appropriate adjustment. RESULTS: We found that the T allele of rs13266634 was protectively associated with lean (OR = 0.810, P = 6.91E-04) but not obese T2D with considerable heterogeneity (P = 0.018). This allele also conferred significant protection with T1D of single (OR = 0.847, P = 9.76E-03), but not multi autoantibodies with substantial heterogeneity (P = 0.005). This variant significantly affected OGTT-related insulin release in lean (P = 2.66E-03, 3.88E-03 for CIR and DI, respectively) but not obese healthy individuals. Furthermore, rs13266634 T allele correlated with the risk of ZnT8A (OR = 1.440, P = 3.31E-05) and IA-2A (OR = 1.219, P = 1.32E-03) positivity, with more effect size in children/adolescents compared with adult-onset T1D subtypes. CONCLUSIONS: These suggested that the SLC30A8 rs13266634 variant might be put into genetic risk scores to assess the risk of the subtypes of T1D and T2D and their related clinical features.

Differences in Maturation Status and Immune Phenotypes of Circulating Helios+ and Helios- Tregs and Their Disrupted Correlations With Monocyte Subsets in Autoantibody-Positive T1D Individuals.

CD4 Tregs are involved in the regulation of various autoimmune diseases but believed to be highly heterogeneous. Studies have indicated that Helios controls a distinct subset of functional Tregs. However, the immunological changes in circulating Helios+ and Helios- Tregs are not fully explored in type 1 diabetes (T1D). Here, we elucidated the differences in maturation status and immune regulatory phenotypes of Helios+ and Helios- Tregs and their correlations with monocyte subsets in T1D individuals. As CD25-/low FOXP3+ Tregs also represent a subset of functional Tregs, we defined Tregs as FOXP3+CD127-/low and examined circulating Helios+ and Helios- Treg subpopulations in 68 autoantibody-positive T1D individuals and 68 age-matched healthy controls. We found that expression of both FOXP3 and CTLA4 diminished in Helios- Tregs, while the proportion of CD25-/low Tregs increased in Helios+ Tregs of T1D individuals. Although the frequencies of neither Helios+ nor Helios- Tregs were affected by investigated T1D genetic risk loci, Helios+ Tregs correlated with age at T1D diagnosis negatively and disease duration positively. Moreover, the negative correlation between central and effector memory proportions of Helios+ Tregs in healthy controls was disrupted in T1D individuals. Finally, regulatory non-classical and intermediate monocytes also decreased in T1D individuals, and positive correlations between these regulatory monocytes and Helios+/Helios- Treg subsets in healthy controls disappeared in T1D individuals. In conclusion, we demonstrated the alternations in maturation status and immune phenotypes in Helios+ and Helios- Treg subsets and revealed the missing association between these Treg subsets and monocyte subsets in T1D individuals, which might point out another option for elucidating T1D mechanisms.

Temporal metabolic and transcriptomic characteristics crossing islets and liver reveal dynamic pathophysiology in diet-induced diabetes.

To investigate the molecular mechanisms underlying islet dysfunction and insulin resistance in diet-induced diabetes, we conducted temporal RNA sequencing of tissues responsible for insulin secretion (islets) and action (liver) every 4 weeks in mice on high-fat (HFD) or chow diet for 24 weeks, linking to longitudinal profile of metabolic characteristics. The diverse responses of α, ß, and δ cells to glucose and palmitate indicated HFD-induced dynamic deterioration of islet function from dysregulation to failure. Insulin resistance developed with variable time course in different tissues. Weighted gene co-expression network analysis and Ingenuity Pathway Analysis implicated islets and liver jointly programmed ß-cell compensatory adaption via cell proliferation at early phase and irreversible islet dysfunction by inappropriate immune response at later stage, and identified interconnected molecules including growth differentiation factor 15. Frequencies of T cell subpopulation showed an early decrement in Tregs followed by increases in Th1 and Th17 cells during progression to diabetes.

The associations between three genome-wide risk variants for serum C-peptide of T1D and autoantibody-positive T1D risk, and clinical characteristics in Chinese population.

AIMS: Recent meta-genome-wide association studies identified several genetic variants associated with beta-cell function in type 1 diabetes (T1D). The aim of this study was to investigate the associations between these variants and T1D risk, C-peptide levels, islet-specific autoantibodies, and lipid levels in Chinese Han population. METHODS: A total of 1005 unrelated autoantibody-positive T1D cases and 1417 healthy controls were included, which were genotyped for rs559047, rs9260151, and rs3135002. T1D individuals were measured for both C-peptide and lipid levels. Logistic regression models were used to examine these associations. RESULTS: We found that rs3135002 A allele showed a genome-wide significant association with T1D risk (OR = 0.22, 95% CI = 0.17-0.30; P = 7.49 × 10-27), and significant heterogeneity of effect size was observed between early-onset and later-onset T1D subgroups (I2 = 80% and P = 0.026). Rs559047 had a nominal association with fasting C-peptide levels in newly diagnosed T1D individuals (P = 0.036). Moreover, rs3135002 A allele was significantly associated with GADA positivity (OR = 0.52, 95% CI = 0.30-0.91, P = 0.02). In addition, nominal correlations were observed with HDL levels for rs559047 (P = 0.042), while LDL levels for rs9260151 (P = 0.032) in T1D individuals. CONCLUSIONS: Our results indicate that there are both similarities and differences for the associations of genetic variants among T1D development, progression, and related autoimmunity, metabolic traits.

High Prevalence of a Monogenic Cause in Han Chinese Diagnosed With Type 1 Diabetes, Partly Driven by Nonsyndromic Recessive WFS1 Mutations.

It is estimated that ∼1% of European ancestry patients clinically diagnosed with type 1 diabetes (T1D) actually have monogenic forms of the disease. Because of the much lower incidence of true T1D in East Asians, we hypothesized that the percentage would be much higher. To test this, we sequenced the exome of 82 Chinese Han patients clinically diagnosed with T1D but negative for three autoantibodies. Analysis focused on established or proposed monogenic diabetes genes. We found credible mutations in 18 of the 82 autoantibody-negative patients (22%). All mutations had consensus pathogenicity support by five algorithms. As in Europeans, the most common gene was HNF1A (MODY3), in 6 of 18 cases. Surprisingly, almost as frequent were diallelic mutations in WFS1, known to cause Wolfram syndrome but also described in nonsyndromic cases. Fasting C-peptide varied widely and was not predictive. Given the 27.4% autoantibody negativity in Chinese and 22% mutation rate, we estimate that ∼6% of Chinese with a clinical T1D diagnosis have monogenic diabetes. Our findings support universal sequencing of autoantibody-negative cases as standard of care in East Asian patients with a clinical T1D diagnosis. Nonsyndromic diabetes with WSF1 mutations is not rare in Chinese. Its response to alternative treatments should be investigated.

Hepatocytes derived extracellular vesicles from high-fat diet induced obese mice modulate genes expression and proliferation of islet ß cells.

Liver secretes proliferative factors participating compensatory hyperplasia of islets during obesity and insulin resistance. Extracellular vesicles (EVs) mediate intercellular communication by delivering inner factors to recipient cells. This study explored the biological effects of hepatocellular EVs on islet ß cells during obesity. Compared with standard chow diet (CD), hepatocellular EVs derived from high-fat diet (HFD) induced obese mice promoted proliferation of ß cell line-MIN6 cells, but didn't influence their insulin secretion. Microarray analysis found 13 miRNAs with significantly differential expression in hepatocellular EVs between HFD with CD group. Meanwhile, RNA-sequencing detected 80 genes with significantly differential expression in MIN6 cells treated with HFD and CD hepatocellular EVs respectively. Six miRNAs and 11 potential target genes were pre-screened by synthesizing TargetScan prediction and RNA-sequencing results. After miRNA mimic transfection and testing the expressions of target genes and cell vitality, miR-7218-5p was verified to affect MIN6 cell proliferation through targeting CD74 gene. SiRNA transfection and dual luciferase reporter assay further confirmed the binding and regulation of miRNA-7218-5p on CD74. These findings suggest HFD induced obesity could change miRNA profiles in hepatocellular EVs, which modulate expression of multiple genes and proliferation of MIN6 cells and maybe mediate compensatory hyperplasia of islets.

Identification of Novel T1D Risk Loci and Their Association With Age and Islet Function at Diagnosis in Autoantibody-Positive T1D Individuals: Based on a Two-Stage Genome-Wide Association Study.

OBJECTIVE: Type 1 diabetes (T1D) is a highly heritable disease with much lower incidence but more adult-onset cases in the Chinese population. Although genome-wide association studies (GWAS) have identified >60 T1D loci in Caucasians, less is known in Asians. RESEARCH DESIGN AND METHODS: We performed the first two-stage GWAS of T1D using 2,596 autoantibody-positive T1D case subjects and 5,082 control subjects in a Chinese Han population and evaluated the associations between the identified T1D risk loci and age and fasting C-peptide levels at T1D diagnosis. RESULTS: We observed a high genetic correlation between children/adolescents and adult T1D case subjects (r g = 0.87), as well as subgroups of autoantibody status (r g ≥ 0.90). We identified four T1D risk loci reaching genome-wide significance in the Chinese Han population, including two novel loci, rs4320356 near BTN3A1 (odds ratio [OR] 1.26, P = 2.70 × 10-8) and rs3802604 in GATA3 (OR 1.24, P = 2.06 × 10-8), and two previously reported loci, rs1770 in MHC (OR 4.28, P = 2.25 × 10-232) and rs705699 in SUOX (OR 1.46, P = 7.48 × 10-20). Further fine mapping in the MHC region revealed five independent variants, including another novel locus, HLA-C position 275 (omnibus P = 9.78 × 10-12), specific to the Chinese population. Based on the identified eight variants, we achieved an area under the curve value of 0.86 (95% CI 0.85-0.88). By building a genetic risk score (GRS) with these variants, we observed that the higher GRS were associated with an earlier age of T1D diagnosis (P = 9.08 × 10-11) and lower fasting C-peptide levels (P = 7.19 × 10-3) in individuals newly diagnosed with T1D. CONCLUSIONS: Our results extend current knowledge on genetic contributions to T1D risk. Further investigations in different populations are needed for genetic heterogeneity and subsequent precision medicine.

Follicular Regulatory T Cells Are Associated With ß-Cell Autoimmunity and the Development of Type 1 Diabetes.

OBJECTIVE: Impaired follicular regulatory T (Tfr) cells enhance T follicular helper cells activity, resulting in the expansion of autoreactive B cells and autoantibody production. However, the role of Tfr cells in the pathogenesis of type 1 diabetes (T1D) is unclear. DESIGN: We evaluated the expression and changes in function of circulating Tfr cells by studying patients with T1D alongside those with type 2 diabetes (T2D), first-degree relatives of T1D patients, and healthy controls. We also investigated the effects of Tfr cells on disease development in nonobese diabetic (NOD) mice and in an adoptive transfer model. RESULTS: Tfr cells were significantly decreased in both patient groups. However, they showed different correlations with fasting C-peptide (C-P) and the area under the curve of blood C-P in patients with T1D and T2D. The frequency of Tfr cells was associated with the number of positive autoantibodies and the titer of glutamic acid decarboxylase autoantibody in T1D patients. Furthermore, Tfr cells decreased significantly after 1 year of follow-up. We also observed Tfr cells in four T1D patients treated with rituximab. After rituximab therapy, the frequency of C-X-C motif chemokine receptor 5 (CXCR5)+ programmed death 1+ Tfr cells was decreased and of CXCR5+ inducible costimulator+ Tfr cells was increased in three patients. We also found that Tfr cells were associated with the development of diabetes in NOD mice and an adoptive transfer model. CONCLUSIONS: Tfr cell deficiency could be involved in the pathogenesis of T1D. Therapy with Tfr cells has potential value for T1D. Modulation of these cells may enhance protective immunity to inhibit autoimmune diabetes.

The association between rs1893217, rs478582 in PTPN2 and T1D risk with different diagnosed age, and related clinical characteristics in Chinese Han population.

OBJECTIVE: To investigate the association between polymorphisms in PTPN2 (rs1893217 and rs478582) and type 1 diabetes (T1D) risk with different diagnosed age, as well as related clinical characteristics in Chinese Han population. METHODS: A total of 2270 Chinese Han individuals (1023 T1D patients and 1247 healthy controls) were genotyped for rs1893217 and rs478582. And 306 newly diagnosed T1D patients were measured for C-peptide levels based on a standard mixed-meal tolerance test. In addition, 40 healthy controls were analyzed for different T cell subsets by multi-color flow cytometry. RESULTS: Neither rs1893217 nor rs478582 showed any association with T1D risk under an additive model. Stratified analysis for T1D diagnosed age revealed that rs1893217, but not rs478582, was significantly associated with T1D patients diagnosed age ≤18 (OR =0.80, 95% CI: 0.67-0.97, p = 0.02). For those diagnosed age >18, neither of them showed any association. We also found that rs1893217 had a higher positive rate of ZnT8A (CC vs. TT carrier, OR = 2.07, 95% CI: 1.07-4.03, p = 0.026) and IA-2A (CT vs. TT carrier, OR = 1.36, 95% CI: 1.02-1.80, p = 0.038). Furthermore, for rs478582, compared with TT, healthy individuals carrying CC/CT carriers had significantly lower frequency and Helios expression of naive Treg subsets (p = 0.049 and 0.048 respectively), but not secreting or activating Treg subsets. In addition, we did not find any association between these two polymorphisms and residual ß-cell function in newly diagnosed T1D patients. CONCLUSIONS: Our results suggest that rs1893217 may increase the risk of early-onset T1D and affect humoral immunity, while rs478582 may affect Treg subsets.

CTLA-4 +49 G/A, a functional T1D risk SNP, affects CTLA-4 level in Treg subsets and IA-2A positivity, but not beta-cell function.

To investigate whether CTLA-4 +49 G/A (rs231775), a tagSNP in Asian, is a functional T1D SNP, we genotyped this SNP with 1035 T1D patients and 2575 controls in Chinese Han population. And 1280 controls measured insulin release and sensitivity based on an oral glucose tolerance test; 283 newly diagnosed T1D patients assayed C-peptide level based on a mixed-meal tolerance test. 31 controls were analyzed for different T cell subsets by multi-color flow cytometry. Under additive model, we found that CTLA-4 +49 G/A was significantly associated with T1D (P = 2.82E-04, OR = 1.25, 95% CI: 1.12-1.41), which was further confirmed by meta-analysis (P = 1.19E-08, OR = 1.65, 95% CI: 1.38-1.96) in Chinese Han population. Although we did not find any association between this SNP and beta-cell function in either healthy individuals or newly diagnosed T1D patients, healthy individuals carrying GG/GA genotypes had lower CTLA-4 expression in naïve or activated CD4 Treg subsets (P = 0.0046 and 0.0317 respectively). A higher positive rate of IA-2A was observed among T1D patients with GG genotype compared with AA (OR = 0.51, 95% CI: 0.30-0.84, p = 0.008). Collectively, CTLA-4 +49 G/A reached a GWAS significant association with T1D risk in Chinese Han population, affects CTLA-4 expression in Treg subsets and subsequently humoral immunity in T1D patients.